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Abstract 


Recent modeling studies have revealed that ice crystal number concentration is one of 
the dominant factors in the effect of clouds on radiation. Since the ice crystal 
enhancement factor and ice nuclei concentration determine the concentration, they are 
both important in quantifying the contribution of increased ice nuclei to global warming. 
In this study, long-term cloud-resolving model (CRM) simulations are compared with 
field observations to estimate the ice crystal enhancement factor in tropical and mid- 
latitudinal clouds, respectively. It is found that the factor in tropical clouds is ~10 3 -10 4 
times larger than that of mid-latitudinal ones, which makes physical sense because 
entrainment and detrainment in the Tropics are much stronger than in middle latitudes. 
The effect of entrainment/detrainment on the enhancement factor, especially in tropical 
clouds, suggests that cloud microphysical parameterizations should be coupled with 
subgrid turbulence parameterizations within CRMs to obtain a more accurate depiction of 
cloud-radiative forcing. 
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1. Introduction 


Clouds affect the global energy cycle via radiation (Chamey 1979; Hartmann el al. 
1992) and therefore provide a platform for cloud microphysics to significantly impact 
climate change (National Research Council 2005; Zeng et al. 2009a). The ice crystal 
enhancement factor, defined as the ratio of ice crystal concentration to ice nuclei (IN) 
concentration, varies from one cloud to another (Pruppacher and Klett 1997). Since the 
factor is important in modeling clouds and radiation (Zeng et al. 2008, 2009b), knowing 
its climatological characteristics (e.g., ensemble average, geographic distribution) is 
imperative for quantifying the contribution of increased IN to global warming. In this 
paper, long-term cloud-resolving model (CRM) simulations are compared with field 
observations to infer the ice crystal enhancement factor at different latitudes. 

a. Ice crystal multiplication 

Ice crystal multiplication in clouds is still perplexing. Ice crystal number 
concentrations often exceed IN concentrations estimated at the cloud-top temperature by 
up-to four orders of magnitude (e.g., Koenig 1963; Mossop et al. 1968; Mossop 1985; 
Hobbs and Rangno 1985, 1990; Blyth and Latham 1993). A riming-splintering 
mechanism, identified by Hallett and Mossop (1974), is one candidate for high ice crystal 
multiplication (Blyth and Latham 1997; Phillips et al. 2001). This mechanism works 
when cloud temperatures are between -3 and -8°C and large droplets (>24 pm diameter) 
as well as relatively fast falling (0.7 m/s) ice particles are present (e.g., Hallett and 
Mossop 1974; Mossop 1985). The mechanism, owing to its strict conditions, cannot be 
used to explain high ice crystal multiplication directly. 
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Blyth and Latham (1997), based on field observations, proposed a multi-thermal 
model of cloud glaciation and revealed that the riming-splintering mechanism contributes 
significantly to ice crystal multiplication via fine cloud dynamical structure. Using a 
three-dimensional (3D) high-resolution cloud model, Ovtchinnikov et al. (2000) 
explicitly simulated the fine dynamic structure and confirmed its importance in ice crystal 
multiplication. 

Due to the effects of fine-scale cloud dynamics on microphysics, a CRM with a 
horizontal resolution of ~1 km must be re-evaluated. Since the model cannot represent 
fine cloud dynamical structure explicitly, it cannot properly simulate ice crystal 
multiplication if no corresponding parameterization is introduced. Hence, it is of interest 
to parameterize the ice crystal multiplication in a CRM. 

b. Entrainment and detrainment in tropical clouds 

Fine cloud dynamical structure has been explored from aircraft observations (e.g., 
Malkus and Scorer 1955; Warner 1970; Paluch 1979; Blyth et al. 1988; Damiani et al. 
2006). Updrafts usually take the form of entraining and detraining thermals, and the 
mixing between a thermal and its surrounding air takes place as a series of discrete events 
rather than continuously (e.g., Austin et al. 1985; Damiani et al. 2006). Since the mixed 
thermals eventually move to the level of zero buoyancy (Raymond and Blyth 1986; 
Taylor and Baker 1991; Emanuel 1994), mixing dominates the fine dynamical structure 
in clouds and therefore affects ice crystal multiplication. 

However, the mixing varies greatly from one geographic region to another. Frequent 
downdrafts, for example, have been observed in the Tropics but not in middle latitudes 
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(e.g., Heymsfield et al. 1978; Wei et al. 1998; Igau et al. 1999). After reviewing the 
aircraft observations from over the past decades, Zipser (2003) concluded that undilute 
updraft cores have not been found in the Tropics but are common in severe storms in 
middle latitudes. Based on this meridional variation in fine cloud dynamic structure, it is 
inferred that the ice crystal enhancement factor in the Tropics is much larger than that in 
middle latitudes (see Sec. 4.a for more discussion). Such variation, in the absence of 
direct observations, can be verified with cloud observations via CRM simulations. 

c. Field observations and CRM simulations 

Recent CRM simulations have revealed that cloud ensembles and radiation are 
sensitive to the ice crystal enhancement factor (Zeng et al. 2008, 2009b). Suppose that a 
CRM is assigned different ice crystal enhancement factors for different simulations. If a 
cloud simulation can be made to match the associated observations, then its assigned 
enhancement factor should be a better inference of the in-situ enhancement factor. This 
diagnostic approach (see Sec. 2 for details) can be applied to estimate the enhancement 
factor over various geographic regions, with the aid of field campaign observations. 

Three field campaigns have provided high-quality cloud observations as well as the 
corresponding large-scale forcing: the Tropical Warm Pool - International Cloud 
Experiment (TWP-ICE), the Kwajalein Experiment (KWAJEX) and the Atmospheric 
Radiation Measurement program’s Spring 2000 Cloud field campaign conducted at the 
Southern Great Plains (ARM-SGP). Comparing the cloud observations from the three 
campaigns with CRM simulations, the ice crystal enhancement factor over those regions 
can be estimated and further analyzed to determine how it varies meridionally. 
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The present study investigates the meridional variation in ice crystal enhancement 
factor. It consists of five sections. In Sec. 2, a CRM is described with special attention on 
how to represent IN and ice crystal enhancement factor. In Sec. 3, CRM simulations over 
three geographic regions are carried out and their results compared with observations to 
infer the ice crystal concentrations in each region. In Sec. 4, a theoretical model is 
introduced connecting the ice crystal enhancement factor and fine cloud dynamic 
structure, which is then used to explain the meridional variation in the enhancement 
factor. Section 5 concludes. 

2. Experiment setup 

A 3D CRM, the Goddard Cumulus Ensemble (GCE) model (Tao and Simpson 1993; 
Tao et al. 2003), is used to simulate clouds and radiation. The model is non-hydrostatic 
and anelastic. It takes account of both absorption and scattering for solar radiation and 
both emission and absorption for infrared radiation. Its cloud-radiation interaction has 
been assessed (Tao et al. 1996). The model parameterizes subgrid-scale (turbulent) 
processes with a scheme based on Klemp and Wilhelmson (1978) and Soong and Ogura 
(1980) and incorporates the effects of both dry and moist processes on the generation of 
subgrid-scale kinetic energy. The model uses a three-class ice formulation based on 
Rutledge and Hobbs (1984) with some modifications (Lang et al. 2007; Zeng et al. 2008) 
for cloud microphysics. It includes the sedimentation of cloud ice (Starr and Cox 1985) to 
better simulate clouds in the upper troposphere. It calculates all scalar variables 
(temperature, water vapor, and all hydrometeors) with a positive definite advection 
scheme (Smolarkiewicz and Grabowski 1990). 
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All of the numerical experiments in this study follow the model setup used in 
previous studies (e.g., Johnson et al. 2002; Xie et al. 2005; Xu et al. 2005; Blossey et al. 
2007; Zeng et al. 2007), which simulated clouds with prescribed large-scale forcing 
derived from field observations. The experiments are 3D, using a 1 km horizontal 
resolution and a vertical resolution that ranges from 42.5 m at the bottom to 1 km at the 
top. The model uses 256 x 256 x 41 gridpoints and a time step of 6 seconds for 
integration. Other model parameters are detailed in Zeng et al. (2008, 2009b). 

The model has five prognostic hydrometeor variables: the mixing ratios of cloud 
water, rainwater, cloud ice, snow and graupel. It parameterizes the Bergeron process via 
the IN concentration and ice crystal enhancement factor. The conversion rate of cloud ice 
to snow due to vapor deposition is expressed as (Zeng et al. 2008, 2009b) 

max[2o,(3^ - m IS( p' x fiN i )m a I J^O], (2.1) 

and the conversion rate of cloud water to ice as 

- + 0 — a 2) m i5oP > ( 2 - 2 ) 

(«2 + 1 )(«2 + 2 ) L J 

where JVJ is the number concentration of active ice nuclei, q t the mixing ratio of cloud ice, 
a\ and <22 the temperature-dependent parameters in the Bergeron process (Koenig 1971), 
p the air density, m/ 5 o=4.8x 10‘ 7 g the mass of an ice crystal 50 pm in diameter, and the 
parameter p represents the ice particle enhancement factor due to riming-splintering and 
other mechanisms (Hallet and Mossop 1974). 

The model uses the formula of Fletcher (1962) to compute the active IN concentration 
in clouds as a function of air temperature T or 

AT, = «,exp[«r 0 -r)] (2.3) 
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where no and (3 are constant. To explore the effects of IN and ice crystal multiplication on 
clouds and radiation, different IN concentrations and ice crystal enhancement factors are 
tested for each large-scale forcing. Table 1 shows the six categories of ice crystal 
concentration used in the numerical simulations. Generally speaking, ice crystal 
concentration increases as the crystal concentration category changes from I to VI. 

Since there were no in-situ observations of IN concentration and ice crystal 
enhancement factor in the three field campaigns, in-situ IN concentrations and crystal 
enhancement factor are implied by comparing CRM simulations with field observations. 
If a CRM simulation with specific values of fj.no and [3 agrees well with field observations 
(i.e., radar, satellite, sounding networks and other measurements), then those values are 
treated as the in-situ ones. 

3. Diagnosing the ice crystal concentrations 

In this section, CRM simulations over three geographic regions are carried out and 
their results compared with cloud observations to diagnose the ice crystal concentrations 
(or the products of IN concentration and ice crystal enhancement factor) in each region. 
The simulations are summarized in Table 2. 

a. TWP-ICE simulations 

To model the tropical clouds observed in TWP-ICE, three simulations T06L, T06M 
and T06H are carried out using the crystal concentration categories II, IV and V (see 
Table 1 for definition), respectively. Obviously, the three simulations correspond to low, 
moderate and high ice crystal concentrations, respectively. Next, their results are 


7 



compared with TWP-ICE observations to determine which category of crystal 
concentration brings about reasonable results. 

TWP-ICE was conducted in Darwin, Australia in January and February of 2006 
during the northern Australian monsoon season (May et al. 2008). It was centered at 
(12°S, 131°E). It provided a great deal of information on clouds such as the liquid and ice 
water contents obtained from ARM Micro-base products (Miller et al. 2003). It also 
provided the large-scale forcing data (e.g., vertical motion and horizontal advective 
tendencies of temperature and moisture) derived using the variational analysis approach 
described in Zhang and Lin (1997) and Zhang et al. (2001). The large-scale forcing data 
represent the mean domain with a center at (12°S, 131°E) and a radius of -120 km (see 
Xie et al. 2009 for more details of the forcing). The present study focuses on the period 
from 2100 UTC 6 to 12 February 2006, a typical monsoon break period, during which 
convection was characterized by intense afternoon thunderstorms with several squall 
lines crossing Darwin in the evening and early morning. 

All three simulations start at 2100 UTC 6 February 2006 and last for 6 days. Figure 1 
displays the three-hour averages of precipitation rate in the simulations. It also displays 
the observed precipitation rate for comparison. Generally speaking, the model captured 
the main precipitation events in spite of some quantitative deviations. The figure also 
shows that the ice crystal concentration can obviously affect precipitation (e.g., day 1.5, 
3,4.5, 5.1 and 5.8). 

Moreover, as the ice crystal concentration increases from low to moderate to high, the 
modeled average precipitation rate decreases while the modeled average precipitable 
water amount increases from 51.5 to 52.5 mm. Compared to the observed precipitable 
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water amount of 54.7 mm, the high ice crystal concentration brings about the most 
reasonable amount of precipitable water, which is consistent with previous results for 
GATE (the Global Atmospheric Research Program’s Atlantic Tropical Experiment) 
simulations (Zeng et al. 2009b). 

In addition, the ice crystal concentration can affect cloud ensembles significantly. 
Figures 2 and 3 display time-pressure cross sections of the simulated ice contents using 
the low and high crystal concentrations, respectively. The two figures show that cloud ice 
content in the upper troposphere increases significantly with increasing ice crystal 
concentration, but snow and graupel do not. 

The modeled ice content can be compared with the observed to determine which 
category of crystal concentration, if any, is close to the in-situ concentration. Figure 4 
displays the vertical profile of six-day mean ice water content retrieved from the radar 
observations over the Darwin station. The retrieval algorithm has been well tested on thin 
non-precipitating clouds but not on thick precipitating clouds (Dong and Mace 2003). 
Figure 4 also displays the vertical profiles of mean modeled ice content (the sum of the 
cloud ice, snow and graupel mixing ratios) for comparison. Obviously, the ice profile 
using the high crystal concentration is close to the retrieved profile above 360 hPa, but 
those using the low and moderate crystal concentrations are not. 

However, all of the modeled ice contents differ greatly from the retrieved below the 
360 hPa level. This difference between the retrieved and modeled contents is not just due 
to model but also the retrieval algorithm. In the retrieval algorithm, it is difficult to 
distinguish ice from liquid water using radar reflectivity, especially in the mixed-phase 
region (X. Dong, personal communication, 2008). Hence, the retrieval in the mixed-phase 
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region may contain a large error and therefore partly explains the difference between the 
retrieved and modeled ice contents below the 360 hPa level. 

Since ice particles, rather than supercooled drops, are common in the upper 
troposphere (e.g., above the 360 hPa level), the retrieval of ice water content is reliable 
there (Dong and Mace 2003). Hence, based on the comparison between the modeled and 
retrieved ice water contents above the 360 hPa level, it is inferred that the ice crystal 
concentration in the Tropics is high. To support this conclusion on ice crystal 
concentration, satellite observations are compared with CRM simulations over another 
tropical region. 

b. KWAJEX simulations 

In contrast to TWP-ICE, KWAJEX was conducted over a tropical open ocean. It was 
centered at (8.8°N, 167.4°E) and took place from July 23 through September 15, 1999. 
Two numerical simulations are carried out to model the clouds observed during 
KWAJEX. Both simulations start at 0600 UTC 24 July 1999 and last for 52 days. The 
model setup is the same as that for the TWP-ICE simulations. 

The two simulations are driven with large-scale forcing derived using the same 
variational analysis approach described in Zhang and Lin (1997) and Zhang el al. (2001), 
and the crystal concentration categories II and V, respectively. The first one, K3DN, uses 
category V and was introduced briefly in Zeng et al. (2008). The second one, K3DL, uses 
category II and its results resemble experiment K3DD described in Zeng et al. (2008). 

K3DN and K3DL correspond to the high and low crystal concentrations, respectively. 
Their modeled precipitation rate, just as shown in Zeng et al. (2008), decreases while 
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precipitable water increases with increasing crystal concentration. The modeled variables 
(e.g., rainfall rate, precipitable water) approach the observed values for the high crystal 
concentration run. Next, satellite observations are compared with the results to show that 
the cloud ensembles from the high crystal concentration simulation are close to the 
observed. 

A precipitation radar (PR) was included on the Tropical Rainfall Measuring Mission 
(TRMM) satellite (Simpson et al. 1988), which flew over the KWAJEX region obtaining 
vertical profiles of radar reflectivity during the campaign. Figure 5 displays the vertical 
profiles of mean and maximum radar reflectivity obtained from the TRMM observations 
over the period. The figure also displays the vertical profiles of simulated radar 
reflectivity sampled when the satellite flew over the campaign site (see Matsui et al. 2009 
for the computational procedure). As shown in the figure, the reflectivity profiles using 
the high crystal concentration are closer to the observed especially in the upper 
troposphere. 

In summary, both the TWP-ICE and KWAJEX simulations show that the modeled 
quantities (i.e., upper-tropospheric cloud ensembles, precipitable water, and precipitation 
rate) are closer to the observed for the high crystal concentration. Simulations over other 
tropical regions [e.g., GATE and SCSMEX (the South China Sea Monsoon Experiment)] 
also support this conclusion that the ice crystal concentration is high in the Tropics (Zeng 
et al. 2008, 2009b). 


c. ARM-SGP simulations 
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Unlike TWP-ICE and KWAJEX, ARM-SGP was centered at (36.6°N, 96.5°W) and 
thus focused on continental clouds in middle latitudes. One ARM-SGP campaign was 
conducted in the spring of 2000. Three simulations, A00H, A00M and A00L, are carried 
out to model the clouds observed during the campaign. All the simulations start at 1730 
UTC 1 March 2000 and last for 20 days. Their model setup follows those for the TWP- 
ICE and KWAJEX experiments except that the surface fluxes are provided by 
observations. 

These experiments with three different ice crystal concentrations (see Table 2 for 
details) have been analyzed previously in Zeng et al. (2009b). Next, the resulting cloud 
ensembles are quantitatively compared with the observed to determine the ice crystal 
concentration during the campaign. 

Figure 6, just like Fig. 4, displays the vertical profile of ice water content retrieved 
from radar observations. The figure also presents the vertical profiles of modeled ice 
water content from the three simulations. As shown in the figure, the ice water content 
from the low ice crystal concentration is closer to the observed above the 400-hPa level, 
which indicates that the ice crystal concentration in spring is rather low in middle 
latitudes. 

Similarly, other ARM-SGP simulations (e.g., those for the summer of 1997 and 2002) 
have shown that the diagnosed ice crystal concentration in summer is quite low in middle 
latitudes too (Zeng et al. 2009a, b). The critical question is why is the climatological 
value of ice crystal concentration for midlatitudinal clouds relatively low but for tropical 
clouds relatively high. This issue is addressed next with theoretical analyses and other 
observations. 
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4. Meridional variation in the ice crystal enhancement factor 

The high ice crystal concentration in tropical clouds is attributed either to a large IN 
concentration or ice crystal enhancement factor. Since IN are usually generated by 
natural and human activities over continents (Pruppacher and Klett 1997; DeMott et al. 
2003; Mohler et al. 2007; Phillips et al. 2008) and there are no IN observations over the 
field campaigns, the in-situ active IN concentrations there are estimated based on their 
geographic locations. 

TWP-ICE, KWAJEX and ARM-SGP occurred in coastal, oceanic and continental 
regions, respectively. Generally speaking, the IN concentrations over the first two 
campaigns should be lower than that over the last. Based on the preceding conclusion that 
the ice crystal concentrations over the tropical campaigns are much higher than that over 
the mid-latitudinal one, it is inferred that the ice crystal enhancement factor in the Tropics 
is much larger than that over middle latitudes, even assuming that the IN concentrations 
over the three campaigns are close to each other. 

The enhancement factor over the tropical campaigns, with the aid of Table 1, is 
quantitatively compared with that over the mid-latitudinal one. The diagnosed crystal 
concentration over ARM-SGP is close to crystal concentration category II, whereas those 
over TWP-ICE and KWAJEX are close to category V. Suppose that the IN 
concentrations (or JVJ) over the three campaigns have the same order of magnitude. Thus, 
the difference in /uNi between category II and V indicates that the enhancement factor // 
over TWP-ICE and KWAJEX is 10 3 -10 4 times larger than that over ARM-SGP, which is 
further explained using a simple theoretical model. 
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a. A theoretical model 


Consider an air parcel with a temperature lower than 0°C and relative humidity of 
100% with respect to water. The parcel fluctuates vertically around its original position 
due to moist turbulence. Let Az m denote the maximum vertical displacement of the parcel 
above its original position. Thus, the ice crystal enhancement factor due to moist 
turbulence is obtained from Eq. (2.3). That is 

fii = exp(fiy£z„) (4.1) 

where y s is the saturated adiabatic lapse rate. Obviously, (4.1) takes no account of the ice 
crystal multiplication due to the riming-splintering mechanism and other multiplication 
processes. 

Expression (4.1) shows that the enhancement factor increases significantly with 
increasing Az m . Suppose J3= 0.6 and ^=7° C/km. Thus, //= 10 when Az m =548 m, and // 
reaches 10 3 when the vertical displacement is 1.5 km or so. 

Convective updrafts and downdrafts, as a part of moist turbulence, can determine the 
maximum vertical displacement of air parcels. Since the drafts vary in frequency from 
one geographic region to another, they can bring about a geographic variation in the ice 
crystal enhancement factor. Aircraft observations show that there are so many downdrafts 
in the Tropics but not in middle latitudes (e.g., Heymsfield et al. 1978; Wei et al. 1998; 
Igau et al. 1999). And, undilute updraft cores have not been found in the Tropics but are 
common in severe storms in middle latitudes (Zipser 2003). Aircraft observations also 
show that ice crystal concentrations are quite small in updrafts cores but large along the 
edges and in downdrafts (e.g., Damiani et al. 2006). All of the observations are 
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consistent with the concept that frequent downdrafts in the Tropics contribute a great deal 
to the high ice crystal concentration (or crystal enhancement factor) there and therefore 
support (4.1). 

b. Effect of the ice crystal enhancement factor 

The present CRM simulations use jiNt to represent ice crystal concentration. Just as 
shown in (2.3) and (4.1), n and vary differently with height. Next, a new TWP-ICE 
simulation is carried out to show whether the enhancement factor and its vertical 
distribution affect cloud ensembles. 

Simulation T06Hu follows the same setup as T06L and T06H except for the ice 
crystal concentration. Suppose that the IN concentration over TWP-ICE is close to that 
over ARM-SGP. Thus, the IN concentration in the simulations is calculated with (2.3) 
using J3= 0.6 and «o=l .2x 10' 9 cm' 3 . 

Based on (4.1) and Blyth and Latham (1997), the enhancement factor is approximated 
with 

H = fJ 0 exp [a(T - T 0 )] when T<T 0 (4.2) 

where <2=0.1 K' 1 and /uo =10 4 are chosen. The enhancement factor in (4.2), which is 
different from the IN concentration, is larger when an air parcel is closer to the freezing 
level. 

T06Hu differs from T06L and T06H only in crystal concentration. Figure 7 displays 
the crystal concentrations versus air temperature used in T06Hu, T06L and T06H. The 
concentration in T06Hu is much higher than that in T06L. It is close to that in T06H in 
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general but increases with height less fast than in H06H. In brief, T06Hu uses an 
enhancement factor of ~10 4 and an IN concentration that is close to the ARM-SGP value. 

Figure 8 displays time-pressure cross sections of the cloud ice, snow and graupel 
mixing ratios in the new simulation. It is compared with Fig. 2 to show whether the 
enhancement factor affects upper-tropospheric ice water content, since T06Hu uses an IN 
concentration similar to T06L but an enhancement factor of ~10 4 . The contrast between 
the two figures indicates that the large enhancement factor brings about a high ice water 
content in the upper troposphere. 

Figure 9 displays the time-pressure cross section of retrieved ice water content over 
the Darwin station. Compared with Figs. 2 and 8, it shows that the modeled cloud ice 
contents in T06Hu as well as in T06H are closer to the retrieved ice water content than 
that from T06L, indicating that the ice crystal enhancement factor over TWP-ICE is ~10 4 . 

Contrasting Figs. 8 and 3 reveals the effect of the vertical distribution of the 
enhancement factor on upper-tropospheric ice water content, because T06Hu and T06H 
use similar crystal concentrations but different vertical distributions of crystal 
concentration. In reference to the retrieved ice water content shown in Fig. 9, the 
modeled cloud ice content in the upper troposphere has a better cloud structure from day 
4 to 6 in T06Hu than in T06H, which implies that the vertical distribution of the 
enhancement factor can significantly affect cloud ensembles and therefore should be 
determined by a physically-based scheme. 

5. Conclusions and discussion 
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Since both the IN concentration and ice crystal enhancement factor strongly impact 
upper-tropospheric ice water content (e.g., Zeng et al. 2008, 2009b), they are both 
important in quantifying the contribution of increased IN to global warming (Zeng et al. 
2009a). Due to the sparse cloud sampling by aircraft, no climatological information on 
the enhancement factor is currently available. In this paper, long-term cloud simulations 
are compared with field observations to infer a meridional variation in the enhancement 
factor, which is summarized as follows. 

• Long-term CRM simulations are compared with TWP-ICE, KWAJEX and ARM- 
SGP cloud observations. It is found that the ice crystal enhancement factor in the 
Tropics is ~10 3 -10 4 times larger than that in middle latitudes. 

• The significant decrease in the enhancement factor with increasing latitude makes 
physical sense. Since fine cloud dynamic structure (e.g., convective downdrafts) 
can affect the enhancement factor greatly and it varies significantly from one 
geographic region to another, the frequent downdrafts or strong vertical mixing in 
the Tropics bring about the large enhancement factor there. 

• Consider CRMs with a horizontal resolution of ~1 km. Since they do not simulate 
moist eddies explicitly, they should parameterize the effect of moist eddies on ice 
crystal concentration. In the present study, CRM simulations with different 
assigned crystal concentrations are compared with observations to infer (or 
diagnose) in-situ crystal concentrations. It is found that the ice crystal 
concentration or ice crystal enhancement factor varies significantly with latitude. 
Hence, to obtain more accurate prediction of cloud radiative forcing in a CRM, 
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the representation of ice crystal concentration should be coupled with subgrid 
turbulence parameterization. This is usually overlooked in current CRMs. 
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Table 1 Crystal concentration categories used in numerical simulations 


Category 

/Mo (cm' 3 ) 

fi 

I 

6.0x1 O' 10 

0.3 

II 

■ 

1 

0.4 

III 


0.45 

IV 


0.5 

V 

1.2x1 O' 6 

0.6 

VI 

1.2x1 O' 5 

0.7 
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Table 2 A list of all the numerical simulations 


Field 

Campaign 

Numerical 

Experiment 

no 

(cm' 3 ) 

P 

F 

Ice Crystal 
Concentration 

Close to 
Observations 


T06H 


0.6 

1.2 

High 

Yes 

TWP-ICE 

T06M 


0.5 

1.2 

Moderate 

No 


T06L 


0.4 

1.2 

Low 

No 


T06Hu 


0.6 

Eq. (4.2) 

High (see text) 

Yes 

KWAJEX 

K3DN 


0.6 

1.2 

High 

Yes 


K3DL 


0.4 

1.2 

Low 

No 


AOOH 


0.6 

1.2 

High 

No 

ARM-SGP 

AOOM 


0.5 

1.2 

Moderate 

No 


AOOL 


0.4 

1.2 

Low 

Yes 
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Caption 

Fig. 1 Time series of surface precipitation rate for the TWP-ICE observations and three 
simulations that start at 2100 UTC 6 February 2006. Thick line represents the 
observations. Dashed, dot-dashed and dotted lines represent the modeling results with 
low (T06L), moderate (T06M) and high ice crystal concentrations (T06H), respectively. 

Fig. 2 Time-pressure cross sections of cloud ice (top), snow (middle) and graupel 
(bottom) mixing ratios for the TWP-ICE simulation T06L with the low ice crystal 
concentration. 

Fig. 3 Same as Fig. 2 except for T06H with the high ice crystal concentration. 

Fig. 4 Six-day mean vertical profiles of ice water content from the TWP-ICE 
observations and the three simulations using the low (T06L), moderate (T06M) and high 
ice crystal concentrations (T06H), respectively. 

Figure 5 Vertical distributions of radar reflectivity from KWAJEX observations and two 
simulations. Blue, red and green symbols denote the observations and simulations K3DL and 
K3DN with the low and high ice crystal concentrations, respectively. Hollow and solid 
symbols represent mean and maximum radar reflectivity, respectively. 
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Fig. 6 Twenty-day mean vertical profiles of ice water content from the ARM-SGP 
observations and the three simulations using the low (AOOL), moderate (AOOM) and high 
ice crystal concentrations (AOOH), respectively. 

Fig. 7 Ice crystal concentration versus air temperature used in simulations T06L (thin 
solid), T06H (dashed) and T06Hu (thick solid line). 

Fig. 8 Same as Fig. 2 except for T06Hu. 

Fig. 9 Time-pressure cross sections of the retrieved ice water content over the Darwin 
station starting at 2100 UTC 6 February 2006. 
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Fig. 1 Time series of surface precipitation rate for the TWP-ICE observations and three 
simulations that start at 2100 UTC 6 February 2006. Thick line represents the 
observations. Dashed, dot-dashed and dotted lines represent the modeling results with 
low (T06L), moderate (T06M) and high ice crystal concentrations (T06H), respectively. 
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Fig. 2 Time-pressure cross sections of cloud ice (top), snow (middle) and graupel 
(bottom) mixing ratios for the TWP-ICE simulation T06L with the low ice crystal 


concentration. 


31 






Cloud Ice (g/kg) 


TWP-ICE (3D) 



Snow (g/kg) TWP-ICE (3D) 


115 r 
-.215 r 
0.315 r 
— 415 r 
2 515 r 
w 615 E- 
a> 715 - 
815 E- 
915 =■ 
1015 


- 1 — i x — i — 1 — i — 1 — i — ■ — r 


1 1 in 





0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.20 0.40 


I ■ I 

4 5 6 

| Time (day) 


Graupel (g/kg) 


TWP-ICE (3D) 


-m 215 

q.315 



0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.20 0.40 


Time (day) 


Fig. 3 Same as Fig. 2 except for T06H with the high ice crystal concentration. 
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TWP-ICE 



Fig. 4 Six-day mean vertical profiles of ice water content from the TWP-ICE 
observations and the three simulations using the low (T06L), moderate (T06M) and high 
ice crystal concentrations (T06H), respectively. 
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Figure 5 Vertical distributions of radar reflectivity from KWAJEX observations and two 
simulations. Blue, red and green symbols denote the observations and simulations K3DL and 
K3DN with the low and high ice crystal concentrations, respectively. Hollow and solid 
symbols represent mean and maximum radar reflectivity, respectively. 
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iwc (g/kg) 


Fig. 6 Twenty-day mean vertical profiles of ice water content from the ARM-SGP 
observations and the three simulations using the low (AOOL), moderate (AOOM) and high 
ice crystal concentrations (AOOH), respectively. 


35 





Fig. 7 Ice crystal concentration versus air temperature used in simulations T06L (thin 
solid), T06H (dashed) and T06Hu (thick solid line). 
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Fig. 8 Same as Fig. 2 except for T06Hu. 
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Fig. 9 Time-pressure cross sections of the retrieved ice water content over the Darwin 


station starting at 2100 UTC 6 February 2006. 
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Popular Summary 

Previous modeling studies have revealed that ice crystal number concentration is one 
of the dominant factors in the effect of clouds on radiation. Since the ice crystal 
enhancement factor and ice nuclei concentration determine the concentration, they are 
both important in quantifying the contribution of increased ice nuclei to global warming. 

In this study, long-term cloud-resolving model (CRM) simulations are compared with 
field observations to estimate the ice crystal enhancement factor in tropical and mid- 
latitudinal clouds, respectively. It is found that the factor in tropical clouds is ~10 3 -10 4 
times larger than that of mid-latitudinal ones, which makes physical sense because 
entrainment and detrainment in the Tropics are much stronger than in middle latitudes. 



The effect of entrainment/detrainment on the enhancement factor, especially in tropical 
clouds, suggests that cloud microphysical parameterizations should be coupled with 
subgrid turbulence parameterizations within CRMs to obtain a more accurate depiction of 
cloud-radiative forcing. 
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